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ABSTRACT 



We use more than 4,500 microflares from the Reuven Ramaty High Energy 
Solar Spectroscopic Imager {RHESST) microflare data set (Christe et ah, 2008, 
Ap. J., 677, 1385) to estimate electron densities and volumetric filling factors of 
microfiare loops using a cooling time analysis. We show that if the filling factor 
is assumed to be unity, the calculated conductive cooling times are much shorter 
than the observed fiare decay times, which in turn are much shorter than the 
calculated radiative cooling times. This is likely unphysical, but the contradic- 
tion can be resolved by assuming the radiative and conductive cooling times are 
comparable, which is valid when the fiare loop temperature is a maximum and 
when external heating can be ignored. We find that resultant radiative and con- 
ductive cooling times are comparable to observed decay times, which has been 
used as an assumption in some previous studies. The inferred electron densities 
have a mean value of 10^^'^ cm~^ and filling factors have a mean of 10~^'^. The 
filling factors are lower and densities are higher than previous estimates for large 
fiares, but are similar to those found for two microfiares by Moore et al. (Ap. J., 
526, 505, 1999). 

Subject headings: Sun: flares — Sun: corona — Sun: activity — conduction — 
radiative transfer 
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INTRODUCTION 



Energy release in flares in solar active regions occurs over many orders of magnitude, 
from large flares to microflares with as low as a millionth the energy content of large flares. 
The latter are A and B GOES class events, which occur more frequently than large flares 



with a negative power law distribution i n number of fla r es as a 



extending over 



Feldman et al. 



many decades iri 



1997 



Wheatland 



2000 



energv (iLin et al. 



Nita et al. 



1984; 



2002 



Dennis 



unction of energy release 



1985 



Paczuski et al. 



Crosbv et al. 



1993 



20051 ) ■ which suggests 



a common energy release mechanism. There are some features in common among flares of 
all si zes, such as radiation in multiple wavelength bands and similar X-ray light curves [see, 
e.g. 



Fletcher et al 



( lin pressi ) for a review]. However, observational differences also exist 



between sm all and large flares. Large flares are asso ciated with higher temperatures than 



small flares ( Feldman et al. 



1995 



Caspi fc Lin 



have steeper spectra than more energetic ones (jChriste et al. 



20101). Also, weaker hard X-ray flares m ay 



2008 



Hannah et al. 



20081 ). 



Statistical studies of hard X-ray microflares have become more comprehensive since 
the laun ch of the Reuve n Ramaty High Energy Solar Spectroscopic Imager {RHESST) 



satellite (iLin et al. 



20021 ). RHESSI SichieYes a lower energy cutoff in the X-ray spectrum 



than previous detectors. By usi ng removable shutter s, RHESSI allows observation of both 



l arge and sma 



( jChriste et al. 



I flar es [see, e.g., 



Hannah et al. 



( I2OIOI ) for a description]. A recent study 



20081) used a new flare-finding technique to identify over 24,000 microflares 



from March 2 002 - March 2007. Statis tical analyses of the microflare properties were 



carried out by 



Christe et al. 



f l2008[ ) and 



Hannah et al. 



( 120081 ). This large data set allows for 



unprecedented studies of flare properties. 

Here, we use this large data set to infer the electron density and the volumetric filling 
factor of the microflare loops in the RHESSI data set. The volumetric filling factor is 
the fraction of the flare loop volume from which radiation is detected. While the filling 
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factor is often thought of as a robust parameter for flare loops, it should be noted that its 
determination is potentially instrument- and re solution-dependent . A previous estimate of 



the density assume d the filling 



analysis [see, e.g., 



Moore et al. 



a ctor w as unity (JHannah et al 



20081 ). We use a cooling-time 



( 1l980l )]. to argue that the volumetric filling factors of 



microfiare loops may be conside rably smaller than u nity, implying densities considerably 



higher than the estimates from iHannah et al. 



(120081 ). The filling factors and densities 



we find are consist ent with a previous study of two microflares observed with Yohkoh 



(JMoore et al 



19991 ) 



The RHESSI microfiare data set is described in Sec. [H The analytical technique is 
reviewed and critiqued in Sec. [31 Results and uncertainty estimates are presented in Sec. HI 
Finally, conclusions are discussed in Sec. [51 



and 



2. OBSERVATIONAL DATA 



A thorough d e script ion of the observational data is given by 



Hannah et al. 



Christe et al. 



(120081) 



(120081 ): the details most salient for the present study are summarized 



here. The data set consists of all microflares observed with RHESSI between March 2002 
and March 2007. The microflares were found to exclusively occur in active regions. The 
events were identified as local maxima in the count rate of 6-12 keV photons having the 
appropriate sign of the time rate of change of the count rate on either side of the maxima 
with signal-to-noise ratio sufficiently large. A total of 24,097 events were identified. Of 
these, spectral fitting and imaging analysis was possible for 4,567 events, which allowed 
for a determination of the plasma parameters required for the present analysis, namely 
the microfiare decay time td, the emission measure EM, the temperature T, and the loop 
length L and width w. 
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Histograms showing distributions of the parameters used in the present study are given 
in Fig. [TJ These and all distributions use 45 equally sized bins. The de cay time T£, of the 



20081 ) . which employs a 



flare is determined from the flare finding algorithm ( IChriste et al.l 
condition on the time derivative of the count rate to determine the end of the event. Their 
distribution is shown in panel (a); the median value is 192 sec. The minimum decay time of 
any microflare in the present study is 4 sec. 

The emission measure EM (assumed isothermal) and temperature T are determined 
from spectral fitting of RHESSI hard X-isiy data to a mo del assuming an isothermal plasma 



with a power law tail above 10 keV (JHannah et al. 



20081 ). The logarithm of EM is plotted 



in (b), with the temperature in (c). Median values of 13.0 MK and 2.38 x 10^^ cm ^. The 



uncertainty in E M was estimated as 



T is less than 1% (jHannah et al. 



20081). 



10% and the estimated (statistical) uncertainty in 



Finally, the physical si z e of t 



described by 



Hannah et al. 



;he flare loops was estimated using visibility forward fitting, 
( 120081 ). which fits several Gaussian sources along the curved loop 
to estimate the full-length L and a central Gaussian FWHM to provide a measure of the 
full- width w of the flare loops. The volume is estimated as V = tt{w/2)'^L. Distributions 
of L,w, and the logarithm of V are plotted in panels (d), (e), and (f), respectively, with 
median values of 2.09 x 10^ cm, 0.665 x 10^ cm, and 7.38 x 10^^ cm^. The estimated 
(statistical) uncertainty in L and w is ~ 20%, which is the standard deviation of repeated 



(100) fit at tempts with t 



le vis ibility amplitudes randomized within their statistical error 



each time ( jHannah et al. 



120081 ). There are also systematic errors such as projection effects 
and the assumption of a circular cross-section of the loops which are not included in the 
estimate. Another possible source of systematic error is that the observations are from 
particles which have considerably higher energy than the thermal background, so the 
determination of L may be an underestimate of the overall size of the loop because of the 



-6- 



absence of thermal particles in the data. 



DATA ANALYSIS 



The loop electron density rig was not measured, but can be estimated using the 
definition of the isothermal emission measure EM^ 

EM = f nldV, (1) 

Jv 

where dV is a differential volume element and V is the total volume of a flare loop. The 

simplest and most common way to estimate Ug is to assume it is uniform over the volume, 

which implies 

(2) 



rip 



V 



This expression is correct if radiation can be detected from all electrons which are present 
in the loop, i.e., the loop is assumed to be optically thin. If flnite optical depth effects are 
present in a particular loop, then the observed X-ray flux from that loop does not include 
direct contributions from all electrons, so the actual density would be higher than the 
estimate in Eq. (|2]). Therefore, Eq. ([2]) provides a lower bound on rig. 



An improvement on this technique comes from deflning the so-called fllling factor 
terms of the fllling factor, the characteristic loop electron density rig is 



Hp 



'EM 



In 



(3) 



Since Eq. (|2]) gives a lower bound on ng, is a positive number between and 1. 

Here, estimates of the density will be tested using a c ooling time analys i s. Simi 



analyses ha.ve been performed previously in many contexts (iMoore et al. 



Stern et al 



1983 



Lin et al 



1992 



Cargill 



1993 



Moore et al 



1999 



1980 



Haisch 



ar 



1983 



Shibata fc Yokovama 



1999 



Aschwanden et al 



200 



Cargill fc Klimchuk 



2004J : 



MuUan et al 



2006 



Jiang et al 
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2006 : iVrsnak et al.l 120061 : iTsiropoula et al.l 120071 : iCassak et al.l 120081 : lAschwanden et al 
20081 ). Cooling time scales are estimated using a scaling analysis (replacing derivatives by 
finite differences of characteristic scales) of the hydrodynamic temperature equation for a 
compressible optically-thin plasma, 

nks dT 



-nkeTV ■ V + kV'T - nf A(T) + Q,,t 



(4) 



7- 1 rft 

where T is the temperature, n is the total plasma density, 7 is the ratio of specific heats, 
V is the bulk flow velocity, ks is Boltzmann's constant, k(T) is the coefficient of thermal 
conductivity, A(T) is the radiative loss function for an optically thin plasma, and Qext is 
the volumetric heating rate from external sources. Here, we assume quasi-neutrality so that 
n ~ 2ne and that the plasma is an ideal gas with 7 = 5/3. Comparing the left hand side to 
the radiative loss term gives a radiative decay time scale t^ of 

3kBT 



-Tr 



neA(T)- 

Comparing the left hand side to the conduction term gives a conductive decay time tc of 

3n,kB{L/2f 



(5) 



Tc 



<T) 



(6) 



where L/2 is half the length of the flare loop which is the distance from loop top to the 
solar surface. 

The radiative loss function is usually taken as a piecewise continuous function 
controlled by different physics at different temperatures. For the temperatures of the 
flare plasmas in the present study [T ~ 8-20 MK from F ig. HTc)], the function al form 
of the radiative loss function is A(r) :^ 5.49 x lO^^V^ Jxhmchuk et al.ll2008l ). For 
the thermal conductivity k(T), we employ t he temperature-depend ent parallel Spitzer 



thermal conductivity oi k(T) = kqT^^'^/ InX ( Spitzer fc HarmI 119531 ). where the coefficient 
kq = 1-84 X 10~^ erg cm^^ s^^ K^^'^, the temperature is in Kelvin, and In A is the Coulomb 
logarithm with A = (3/2e^)(/c^T^/7™e)^/^ for a pure hydrogen plasma. 



This type of scaling analysis has been used previously to estimate cooling times of 
flare loops and determine which mechanism dominates the cooling, typically for large flares. 



Early studies ( JAntiochos &: Sturrock 



19761 ) suggested condu ctive cooling is more efficien t 



but the effects of chromospheric e vaporation slow it down (JAntiochos fc Sturrock 



The role of radiation was studied (JAntiochos 



19781). 



19801), and for a while it was belie ved that 



radiation and conduction act comparably to cool flare loops (jMoore et al. 



19801) because 



the predicted times from the scaling analysis were comparable to observed flare loop decay 
times T£,. From the theoretical perspective, it is reasonable that these time scal es are 
comparable due to the 



1980 



Moore et al. 



1980h . 



unction of the chromosphere as a reservoir for the corona (^Sturrock 



The observational and theoretical result that tq ^ t^ r^ t^j prompted authors to assume 



this relation to esti mate parameters for stef 



were not available (JHaisch 



1983 



Stern et al 



ar fla res for which well-resolved optical data 



19831 ). A recent study comparing predictions 



using this r aodel to independe ntly derived parameters (n, L, T) of stellar loops found good 



agreement (JMuUan et al. 



20061 ) ■ lending credence to the validity of this assumption. 



However, caution must be used in interpreting the scaling analysis time scales as 



genuinely representative of the decay of flare loops. iDoschek et al. 



(119821 ) used simulations 
to suggest that conduction dominates early in time when the temper ature is highest, 



Cardlll (11993h used 



followed by comparable contributions from radiation and conduction. 
a model in which strictly conductive cooling occurred at early times before transitioning 
to radiative cooling at flare maximum because of chromospheric evaporation enhancing 
radiative cooling at late times. Therefore, there need not be a single dominant mechanism 
throughout the duration of the event. 

Another important issue is that the parameters that go into Eqs. ([5]) and ([6]) are tacitly 
assumed to be constant and uniform, but the temperature changes as the loop cools and 
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thus the scaling results are not applicable to finding the time it takes to cool from one 



temperature to another (ICargill et al. 



woul d require time integration (jCulhane et al. 



1995L 



aking into account the change in temperature 



1970 : 



Svestka 



1987 : 



Aschwanden &: Tsiklauri 



20091 ). Since the difference between actual cooling times and the scaling result can be 



significant, the sea. 



( ICargill et al. 



19951 ). 



ing analysis time scales only indicate instantaneous time scales of cooling 



The scaling analysis is on firmer theoretical footing at peak flare temperature. At peak 
temperature, the loo p goes from being h eated to cooling, so the left hand side of Eq. (jl]) is 



zero instantaneously ([Aschwanden 



20071 ). Ignoring the expansion term (which is safe when 



flow speeds are subsonic) and assuming that there is no external heating, the right-hand 
side implies that conduction balances radiation instantaneously at the temperature peak, 
i.e., 

Tc ~ rn. (7) 

This is the same relation as before, but with a very different interpretation. For the 
purposes of the present study, we subscribe to the latter interpretation of expecting equality 
only at peak flare temperature rather than interpreting the scaling results as predictions 
for the actual decay time. (However, a relation to the decay time will be discussed in the 
following section.) Thus, the cooling times are evaluated at the beginning of the cooling 
process. 

It is important to point out aspects left out of the present model that have been 
discus sed in previous s tudies. The scaling analysis ignores p ressure variation along the 



tube (jSerio et al 



198ll ). rad iative cooling at lo op footpoints (lAntiochos fc Sturrock 



1982), 



chromospheric eyapor ation (ICargill et al.lll995f 



Forbes fc Acton 



19961). spatial nonunifo rmity 



multiple loops ( jReeves fc Warren 



, shrinkage of loops 



flAntiochos et al. 



jSvestka et al.i 119871 : 



2000h. and the e ffect of 



2OO2I ) . See lAschwanden fc AschwandenI (120081 ) for an 
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approach incorporating fractal dimensi onal filling of flare loops. A recent study emphasizes 



the role of enthalpy flux in flare loops (JBradshaw fc Cargill 



20101) • If any of these aspects of 



solar flare evolution play a significant role in determining the time-scales of flare decay, the 
results of the present paper could change in detail. 

4. RESULTS 



The density estimates presented in iHannah et al.l ( 120081 ) employed Eq. ([2]) which 



assumes a filling factor of unity. The distribution of the logarithm of Ue under this 
assumption for the data in the present study is plotted in Fig. [2](a). The mean electron 
density of the distribution in Fig. [2|^a) is rig — 10^'* cm~^. As noted earlier, this mean 
density is a lower bound on the true mean density of the flares in the present study. 

Using the density derived for each individual event, we estimate the radiative and 
conductive cooling times using Eqs. ^ and (jH]) for each microflare. The raw values for 
the distributions of the logarithms of the resultant radiative and conductive cooling times 
are plotted in Figs. [2]^b) and (c), respectively. The median values of the calculated tr and 
Tc are 2.06 x 10^ sec and 5.43 sec, respectively. Figs. [2](d) and (e) show the same values 
normalized to the flare's observed decay times r^. Panels (f) and (g) show a scatter plot 
of the logarithm of the calculated cooling times compared to the logarithm of their decay 
times, showing that the time scales are essentially uncorrelated. Panel (d) shows that the 
radiative cooling times are distributed around a peak about 100 times longer than td, while 
panel (e) shows that the conductive cooling times are distributed around a peak almost 100 
times smaller than To-, i.e., tc '^ td '^ tr. This strongly contradicts the hypothesis in 
Eq. ([7]). Note, this assessment assumes that it is reasonable to compare the measured decay 
time to instantaneous e-folding times from the scaling. Since they may differ, this could 
introduce systematic errors in the comparison. However, we suspect that it is not enough 
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to account for the large separation in scales inferred here. 

Assuming the plasma parameters we use (including the density) are correct, it is 
difficult to envision a physical explanation which could justify the disparity in time scales 
because one expects a flare to decay on a time-scale determined by the shortest available 
dissipation time. From this perspective, it is difficult to understand how, in the presence of 
strong conductive cooling, the actual decay time can be much longer than the conductive 
time-scale. A likely cause is the assumption that = 1, which we now relax. 

In order to address the four orders of magnitude disparity which appears to exist 
between the conductive and radiative decay times, we assume Eq. ([7]) holds and use it to 
solve for the density and filling factor. We then check to see if this assumption helps us 
arrive at an internally consistent set of flare decay time scales. Equating Eqs. (^ and (E]) 

and solving for Uf. gives 

Uk(T)T 
n^ — ■ 



This is equivalent to the classical analysis of 



'-MTV 



Rosner et al 



(Il978[ ). Since k is a (weak) 



function of density due to the Coulomb logarithm, we employ an iterative technique to 
self-consistently solve for the density. The procedure is to assume = 1 to obtain a zeroth 
order estimate n^ which is used to calculate the zeroth order kq. The next order of density 
Ui is then determined from the previous kq. This is continued until convergence. We flnd ng 
is determined to eight signiflcant flgures after ten iterations and that the iterative procedure 
changes A(T) by only 10%. Using this value of the density, the flUing factor is obtained 
from Eq. 02]) and cooling times are obtained from Eqs. ^ and (jH]). 

The results of this analysis are displayed in Fig. [31 Panel (a) shows the distribution 
of the logarithm of the radiative cooling time r^, while (b) shows the distribution for the 
logarithm of tr normalized to the flare decay time To- The conductive cooling times tc are 
equal to tr by construction. The median value of tr is 325 sec, which is within a factor of 
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1.7 of the median value of T£,. This is reiterated in panel (c), which is a scatter plot of tr 
and T£). The values do not appear to be correlated, but they are clearly of the same order. 
It is surprising that the scales of the cooling times are essentially equal to a key empirical 
time scale, td'-, nothing in the model requires that such a similarity should emerge from 
the analysis. Of course, there are uncertainties associated with the comparison of observed 
decay times and predicted e-folding scaling times, but the present results lend observational 
support that a model with tc ^ tr r^ To (as has been often assumed before) is consistent, 
at least for the present data set. 

Panel (d) shows the logarithm of the resultant values for the calculated densities from 
Eq. ([8]). The distribution has a mean of rig ~ 10^^'^ cm~^, which is nearly two orders of 



magnitude higher than the reported values in 



Hannah et al. 



(120081 ) and plotted in Fig.EJ^a). 



The inferred filling factors, the logarithm of which is shown in panel (e), have a mean of 
(p ~ 10~^'^. We discuss these results further in the following section. 

Using the calculated densities, one can calculate other properties of the flare loops. 
The logarithm of the total thermal energy Wt = SrieksTeV in the flare loops is shown 
in Fig. m^a), with a median value of 1.57 x 10^° erg. The distribution of the logarithm 
of the calculated gas pressures P = 2nekBTe is shown in Fig. Hl^b), with a median value 
of 1.43 X 10^ erg cm~^. Figure Hl^c) is a scatter plot of the logarithm of the calculated 
fllling factor (p against the electron temperature Tg. The flUing factor is smaller for higher 
temperature loops, which is a consequence of Eqs. ([3]) and ([8]). 

We now discuss the effect of uncertainties in the present analysis. As discussed in 



Sec. [21 the statistical errors in the 



approximately 20%, 10%, and 1% (IChriste et al. 



ength L, emission measure EM, and te mperature T are 



2008; 



Hannah et al. 



20081 ). Standard error 



propagation techniques imply uncertainties for calculated cooling times of 40%, densities 
20%, and filling factors 40%, which are sizable but not unreasonably large. 
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A few potential sources of systematic errors have been noted. They include assuming 
that the measured decay time tj:, corresponds to an e-folding decay time from a scaling 
an alysis, that the repre sentative loop length is assumed to be equal to the values obtained 



by 



Hannah et al. 



(120081 ) determined from particles at energies far above thermal energies, 



and that there is no external heating at the time of peak flare temperature. 

5. DISCUSSION 

In this paper, RHESSI microfiare data are used to estimate the volumetric filling 
factor and the electron density of microfiare loops using an analysis of cooling times. If the 
filling factor is assumed to be unity, then the conductive cooling time of the loop is much 
smaller than the observed decay time, which itself is much smaller than the radiative decay 
time. This is difficult to justify physically. Alternately, if one invokes the hypothesis that 
the radiative and conductive cooling times are comparable at the moment when the flare 
temperature passes through its maximum value (and that cooling due to expansion and 
flare heating are negligible at that time), one can solve for the filling factor and density. 
Mean values for the whole distribution are ~ 10^^'^ and rie ~ 10^^'^ cm~^. Our weakest 
assumption is that fiare heating stops at the peak time of hard X-rays. Since the hard 
X-ray time profile is a convolution of heating and cooling, heating does not necessarily stop 
at the hard X-ray peak time. If heating is present during the decay, even at a low level, 
the cooling times could be longer than derived for the case without heating, and the filling 
factor could be larger than derived here. 



Our estimate of mean densities are higher than those reported by 



Hannah et al 



(120081). 



The authors are aware of only one other systematic study of microfiare loop densities 



or filling factors ( iMoore et al. 



19991 ). who used Yohkoh to study two microfiare strands. 



Using an identical analytical technique as the present study, they found densities between 
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Hf. ~ 



10^° cm ^ and 10^^'^ cm ^, with filling; factors between 



the values in the two microfiares studied by (JMoore et al 



0-3-2 ^^^ ^Q-2.8_ rp^^ 



19991 ) from Yohkoh are in good 



agreement with results from the large RHESSI data set studied here. 



We now compare the present results with observations of l arge flares. ICul 



(jl994j ) found ~ 1 and rig ~ 3 x 10^^ cm ^ for an M-class flare 
(/)^ .01-0.2 and n, - 7 x 10^-1.5x10^° cm'^ for a C-class flare, 
(I2OOII ) found ~ 1 and ng ~ 1.5 x 10^° cm'^ for the Bastille Day flare, 



Varadv et al. 



lane et al. 



(2000) found 



Aschwanden fc Alexander 



Teriaca et al. 



(120061 ) 



found ~ 0.2 — 0.5 and rie ~ 10^*^ cm ^ for a C-class flare, and iRaymond et al.i (i2007i ) found 
(j) > 0.01 and n^ ~ 10^^ c m"^ for X-class fla res. Other examples of fillin g factors include 



20091 ). 0.0 4-0.07 in coronal ho 



0.3 in active region loops ( jLandi et al. 

20091 ). and near unity in many coronal hole jets (iDoschek et al. 



study (IChifor et al 



es (Abramenko et al. 



2OIOI) but 0.03 in another 



20081 ) . As noted earlier, the determination of filling factors can depend 
on detector resolution and wavelength. Nonetheless, the filling factors obtained here for 
microfiares are at least 10 times smaller than those reported for large flares, which is likely 
statistically significant. Densities are slightly higher for microfiares in the present study 
than for larger flares in previous studies. 

Given the characteristics of the microflare data set considered here compared to large 
flares, it is perhaps not surprising that the filling factors are small. The microflare loops 
in the present study occur in active regions, just as large flares do. The mean sizes of 
the loops in the present study are comparable to those of larger flares. The loops have 
energies of rieksTV ~ 10^*^ ergs deposited into them by the flare (using characteristic values 
from Fig. [1]), which is about 10^ times smaller than large flares. Thus, the loops are of 
similar size but acquire less energy, which could lead to a smaller filling factor. To estimate 
the size of the region for which radiation is detected, we note the radiating volume is 
K = (pV. Assuming that the length of the radiating plasma is L, the effective width w^, of 
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the radiating plasma is given by V^, ~ 7r(u;^,/2)^L. For the present parameters, this imphes 
loops with total thickness w^ ~ (p^/'^w ~ u;/100 ~ 4 x 10^ cm, which would imply there is 
an unseen substructure of thin strands within the flare loops. 



Various 



e.g. 



ines of evidence indicate that there are smaller-scale structures in the corona, 
MuUanl ( 1l990l ). Radio polarizat ion data point to the existence of structures in the 



corona which are ~ 100 km in size (JMelrose 



19751 ). The possibility that 100 km structures 



are associated with collapsing magnetic reconnection sites in the corona was discussed by 



MuUanl (jl980l ): using constraints on the collapse time-scales and coronal Alfven speeds. 



transverse dimensions of order 100 km were found to be typical of reconnection sites in 
the corona. More recently, there is evidence that X-class flare loops are composed of thin 



threads from high resoluti on observations, with structure at scales o : 



(1" :ii 10*^ cm) and below ( iDennis &: Pernak 



2009 



Kontar et al. 



2010 : 



a few arcseconds 



Krucker et al. 



2010l ) 



There is also abundant evidence from the footpoints of flaring loops that most of t he 



emission i s spat ially unresolved, such as in TRACE white-light flares (JHudson et al. 



Xu et al. 



20061). 



(120061 ) found a core region within a halo region in two X-class white-light flares. 



reporting a ratio of the area of the core to the halo of 4% and 25%, respectively. Also, 
simulations of loops comprised of many small scale filamen ts were able to repr oduce cooling 



characteristics of large flare loops observed with TRACE (IWarren et al. 



20031). Thus, the 



conclusion that there is small substructure of flare loops is not without precedent. 



The pr ediction of sma l 
flare loops. iHannah et al. 



scale loops has implications for the heating mechanism of the 



( 120081 ) estimated that the non-thermal power in accelerated 
electrons during the time of peak emission in the RHESSI microHaies is 10^^ erg s~^. 
For loops of area 10^^ cm^ as predicted by the present results, the energy deposition rate 
per unit area would be 10^^ erg s^^ cm"^. This is an enormous value, as discussed in 



Krucker et al. 



( 120101 ). Hence, if the fllling factor is indeed ~ 10 ^, then microflares are 
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unlikely heate d by electron beams. A r ecent model that the flare energy is transported by 



Alfven waves (JFletcher fc Hudson! |2008| ) would not be ruled out by the data 



An interesting result of the present study is that the conductive and radiative cooling 
times derived by assuming their equality at the time of maximum temperature are 
comparable to the observed microflare decay times. A possible ramiflcation of this result 
is that it lends credence to the assumption that the conductive and radiative times are 
comparable to the decay time, tr r-^ tq ^ to- This is relevant to stellar flare studies in 
which plasma parameters were obtained under such a ssumptions 



aisch 



19881 : 



Stern et al. 



20061 ) found results of 



19831 ). In addition, a previous study of stellar flares ( iMuUan et al. 
this model are largely consistent with independently determined plasma parameters. If one 
believes that the scaling analysis cooling times actually represent physical cooling times for 
loops, the result suggests that conductive and radiative cooling act at comparable levels to 
cool flare loops, at least for the microflares in the present study. 

Future work could include efforts to incorporate physical effects left out of the model 
as summarized at the end of Sec. O Also, future studies could further try to minimize the 
systematic errors discussed in Sec. HI These can be addressed both with observations and 
with numerical modeling. Also, the study of the cooling times of individual events will help 
determine the validity of the cooling time analysis. 

The authors thank G. Holman and B. Dennis for helpful conversations. RNB and PAC 
gratefully acknowledge support by NSF grant PHY-0902479, NASA's EPSCoR Research 
Infrastructure Development Program and the West Virginia University Faculty Senate 
Research Grant program. DJM is supported in part by the Delaware Space Grant. IGH 
is supported by a STFC rolling grant and by the European Commission through the 
SOL AIRE Network (MTRN-CT-2006-035484). RPL is supported in part by the WCU grant 
(No. R31-10016) funded by the Korean Ministry of Education, Science and Technology. 



-17- 



HSH, RPL, and SK are supported through NASA contract NAS 5-98033 for RHESSI. 



-18- 

REFERENCES 

Abramenko, V., Yurchyshyn, V., & Watanabe, H. 2009, Solar Phys., 260, 43 

Antiochos, S. K. 1980, Ap. J., 241, 385 

Antiochos, S. K., DeLuca, E. E., Golub, L., & McMullen, R. A. 2000, Ap. J., 542, L151 

Antiochos, S. K. & Sturrock, R A. 1976, Solar Phys., 49, 359 

— . 1978, Ap. J., 220, 1137 

— . 1982, Ap. J., 254, 343 

Aschwanden, M. J. 2007, Adv. Space Res., 39, 1867 

Aschwanden, M. J. & Alexander, D. 2001, Solar Phys., 204, 93 

Aschwanden, M. J. & Aschwanden, P. D. 2008, Ap. J., 674, 544 

Aschwanden, M. J., Stern, R. A., & Gudel, M. 2008, Ap. J., 672, 659 

Aschwanden, M. J., Tarbell, T. D., Nightingale, R. W., Schrijver, C J., Title, A., 
Kankelborg, C. C, Martens, P., & Warren, H. P. 2000, Ap. J., 535, 1047 

Aschwanden, M. J. & Tsiklauri, D. 2009, Ap. J. Supp. Ser., 185, 171 

Bradshaw, S. J. & Cargill, P. J. 2010, Ap. J., 710, L39 

Cargill, P. J. 1993, Sol. Phys., 147, 263 

Cargill, P. J. & Klimchuk, J. A. 2004, Ap. J., 605, 911 

Cargill, P. J., Mariska, J. T., & Antiochos, S. K. 1995, Ap. J., 439, 1034 

Caspi, A. & Lin, R. P. 2010, Ap. J. Lett., submitted 



-19- 

Cassak, P. A., Mullan, D. J., & Shay, M. A. 2008, Ap. J. Lett., 676, L69 

Chifor, C, Young, P. R., Isobe, H., Mason, H. E., Tripathi, D., Hara, H., & Yokoyama, T. 
2008, Astron. Astrophys., 481, L57 

Christe, S., Hannah, I. G., Krucker, S., McTiernan, J., & Lin, R. P. 2008, Ap. J., 677, 1385 

Crosby, N. B., Aschwanden, M. J., & Dennis, B. R. 1993, Solar Phys., 143, 275 

Culhane, J. L., Philhps, A. T., Indakoide, M., Kosugi, T., Fludra, A., Kurokawa, H., 

Makishima, K., Pike, C D., Sakao, T., Sakurai, T., Doschek, G. A., & Bentley, R. D. 
1994, Solar Phys., 153, 307 

Culhane, J. L., Vesecky, J. P., & Phillips, K. J. H. 1970, Solar Phys., 15, 394 

Dennis, B. R. 1985, Solar Phys., 100, 465 

Dennis, B. R. & Pernak, R. L. 2009, Ap. J., 698, 2131 

Doschek, G. A., Boris, J. P., Cheng, C.-C, Mariska, J. T., & Oran, E. S. 1982, Ap. J., 258, 
373 

Doschek, G. A., Landi, E., Warren, H. P., & Harra, L. K. 2010, Ap. J., 710, 1806 

Feldman, U., Doschek, G. A., & Klimchuk, J. A. 1997, Ap. J., 474, 511 

Feldman, U., Laming, J. M., & Doschek, G. A. 1995, Ap. J., 451, L79 

Fletcher, L., Dennis, B. R., Hudson, H. S., Krucker, S., Phillips, K., Veronig, A., Battaglia, 
M., Bone, L., Chen, Q., Gallagher, P., Grigis, P. T., Ji, H., Liu, W., Milligan, R. O., 
& Temmer, M. 2011, Space Sci. Rev., in press 

Fletcher, L. & Hudson, H. S. 2008, Ap. J., 675, 1645 

Forbes, T. G. & Acton, L. W. 1996, Ap. J., 459, 330 



-20- 

Haisch, B. M. 1983, in lAU Colloq. 71, Activity in Red-Dwarf Stars, ed. P. B. Byrne & 
M. Rodono (Dordrecht: Reidel), 255 

Hannah, I. G., Christe, S., Krucker, S., Hurford, G. J., Hudson, H. S., & Lin, R. P. 2008, 
Ap. J., 677, 704 

Hannah, I. G., Hudson, H. S., Battagha, M., Christe, S., Kasparova, J., Kundu, M., 
Krucker, S., & Veronig, A. 2011, Space Sci. Rev., in press 

Hudson, H. S., Wolfson, C. J., & Metcalf, T. R. 2006, Solar Phys., 234, 79 

Jiang, Y. W., Liu, S., Liu, W., & Petrosian, V. 2006, Ap. J., 638, 1140 

Khmchuk, J. A., Patsourakos, S., & Cargill, P. J. 2008, The Astrophysical Journal, 682, 
1351 

Kontar, E. P., Hannah, I. G., Jeffrey, N. L. S., & Battagha, M. 2010, Ap. J., 717, 250 

Krucker, S., Hudson, H. S., Benz, A. O., Csillaghy, A., & Lin, R. P. 2010, Ap. J., submitted 

Landi, E., Miralles, M. P., Curdt, W., & Hara, H. 2009, Ap. J., 695, 221 

Lin, J., Zhang, Z., Wang, Z., & Smartt, R. N. 1992, Astron. Astrophys., 253, 557 

Lin, R. P., Dennis, B. R., Hurford, G. J., Smith, D. M., Zehnder, A., Harvey, P. R., Curtis, 
D. W., Pankow, D., Turin, P., Bester, M., Csillaghy, A., Lewis, M., Madden, N., van 
Beek, H. P., Appleby, M., Raudorf, T., McTiernan, J., Ramaty, R., Schmahl, E., 
Schwartz, R., Krucker, S., Abiad, R., Quinn, T., Berg, P., Hashii, M., Sterling, R., 
Jackson, R., Pratt, R., Campbell, R. D., Malone, D., Landis, D., Barrington-Leigh, 
C. P., Slassi-Sennou, S., Cork, C, Clark, D., Amato, D., Orwig, L., Boyle, R., Banks, 
I. S., Shirey, K., Tolbert, A. K., Zarro, D., Snow, P., Thomsen, K., Henneck, R., 
Mchedlishvili, A., Ming, P., Fivian, M., Jordan, J., Wanner, R., Crubb, J., Preble, J., 



-21 - 

Matranga, M., Benz, A., Hudson, H., Canfield, R. C, Holman, G. D., Crannell, C, 
Kosugi, T., Emslie, A. G., Vilmer, N., Brown, J. C, Johns-Krull, C, Aschwanden, 
M., Metcalf, T., & Conway, A. 2002, Solar Phys., 210, 3 

Lin, R. R, Schwartz, R. A., Kane, S. R., Felling, R. M., & Hurley, K. C. 1984, Ap. J., 283, 
421 

Melrose, D. B. 1975, Solar Phys., 43, 79 

Moore, R., McKenzie, D. L., Svestka, Z., Widing, K. G., Antiochos, S. K., Dere, K. P., 

Dodson-Prince, H. W., Hiei, E., Krall, K. R., Krieger, A. S., Mason, H. E., Petrasso, 
R. D., Pneuman, G. W., Silk, J. K., Vorpahl, J. A., & Withbroe, G. L. Solar Flares 
(Colorado Associated University Press) 

Moore, R. L., Falconer, D. A., Porter, J. G., & Suess, S. T. 1999, Ap. J., 526, 505 

Mullan, D. J. 1980, Ap. J., 237, 244 

— . 1990, Astron. Astrophys., 232, 520 

Mullan, D. J., Mathioudakis, M., Bloomfield, D. S., & Christian, D. J. 2006, 
Ap. J. Supp. Ser., 164, 173 

Nita, G. M., Gary, D. E., Lanzerotti, L. J., & Thomson, D. J. 2002, Ap. J., 570, 423 

Paczuski, M., Boettcher, S., & Baiesi, M. 2005, Phys. Rev. Lett., 95, 181102 

Raymond, J. C, Holman, G., Ciaravella, A., Panasyuk, A., Ko, Y.-K., & Kohl, J. 2007, 
Ap. J., 659, 750 

Reeves, K. K. & Warren, H. P. 2002, Ap. J., 578, 590 

Rosner, R., Tucker, W. H., & Vaiana, G. S. 1978, Ap. J., 220, 643 



-22- 

Serio, S., Peres, G., Vaiana, G. S., Golub, L., & Rosner, R. 1981, Ap. J., 243, 288 

Shibata, K. & Yokoyama, T. 1999, Ap. J., 526, L49 

Spitzer, L. & Harm, R. 1953, Phys. Rev., 89, 977 

Stern, R. A., Underwood, J. H., & Antiochos, S. K. 1983, Ap. J., 264, L55 

Sturrock, P. A. 1980, Solar flares: A monograph from SKYLAB Solar Workshop II 
(Colorado Associated University Press) 

Svestka, Z. 1987, Solar Phys., 108, 411 

Svestka, Z. F., Fontenla, J. M., Machado, M. E., Martin, S. F., & Neidig, D. F. 1987, Solar 
Phys., 108, 237 

Teriaca, L., Falchi, A., Falciani, R., Cauzzi, G., & Maltagliati, L. 2006, Astron. Astrophys., 
455, 1123 

Tsiropoula, G., Tziotziou, K., Wiegelmann, T., Zachariadis, T., Gontikakis, C., & Dara, H. 
2007, Solar Phys., 240, 37 

Varady, M., Fludra, A., & Heinzel, P. 2000, Astron. Astrophys., 355, 769 

Vrsnak, B., Temmer, M., Veronig, A., & Karhcky, M. 2006, Solar Phys., 234, 273 

Warren, H. P., Winebarger, A. R., & Mariska, J. T. 2003, Ap. J., 593, 1174 

Wheatland, M. S. 2000, Ap. J., 536, L109 

Xu, Y., Cao, W., Liu, C., Yang, G., Jing, J., Denker, C., Emslie, A. G., & Wang, H. 2006, 
Ap. J., 641, 1210 



This manuscript was prepared with the A AS lATfrjX macros v5.2. 



-23- 



(a) 



0.20 



>.0.15 



c 



0.10 

0.05 
0.00 



i 








+ 








+ 


+ 








++ 
++ 


-+++++ 


++-I-^W 



(b) 



0.05 
0.04 

0.03 

0.02 

0.01 
0.00 






-^- 



--#-- 









++ 



+ 



(c) 0.08 

^ 0.06 

D- 0.04 

>^ 0.02 
0.00 







^++ 










+ 
+ 


+ 
-i 










+ 


+ 








+ 




+ 
+ 
+ 










+ 
:++ 



(d) 



o 

c 

a- 

03 



300 600 900 1200 
Decay Time [sec] 

(e) 



44 45 46 47 48 



8 10 12 14 16 18 20 



loQio (Emission Measure [cm ]) Flare Temperature [MK] 



U.IU 
0.08 

0.06 

0.04 

0.02 
0.00 


^ 








+ :+ 
+ ; 








; + 








: + 

+ : ^ 




H*MHf. 


Hutiu4 



0.12 
0.10 
0.08 
0.06 
0.04 

0.02 
0.00 



++t, 



(f) 



o 

c 



0.06 
0.05 
0.04 
0.03 
0.02 

0.01 
0.00 



^<^ 



14 



— ^* — 
+ 



+ 



* 
^"" 






2 4 6 8 10 
Loop Length [1 0^ cm] 







1 



Loop Width [10^ cm] 



24 25 26 27 28 29 
logio(Loop Volume [cm^]) 



Fi g. 1. — Distributions of ra w data used in the p resent study using the methods described 



by 



Christe et al 



(120081 ) and 



Hannah et aL 



(120081 ). Plotted are (a) flare decay time td, (b) 
the logarithm of the emission measure EM, (c) the logarithm of the temperature T, (d) flare 
loop full-length L, (e) flare loop full-width w, and (f ) the logarithm of the flare loop volume 
V. 
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Fig. 2. — Flare parameters assuming a filling factor of unity. Plotted are distributions of 
the logarithms of (a) electron density rig, (b) calculated values of the resistive cooling time 
tr and (c) the conductive cooling time tq, (d) the ratio of radiative to flare decay times 
tr/td) and (e) the ratio of conductive cooling times to flare decay times tc/td. Panels (f) 
and (g) are scatter plots of the logarithms of tr vs. td and tq vs. r^, respectively. 
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Fig. 3. — Flare parameters assuming radiative and conductive cooling times are equal. 
Plotted are distributions of the logarithms of (a) the radiative cooling time tr and (b) the 
ratio of radiative to flare decay times tji/t£,. The conductive cooling time is equivalent to 
the radiative cooling time by construction. Panel (c) shows a scatter plot of the logarithms 
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Fig. 4. — Distributions of the logarithms of (a) the calculated thermal energy Wt = 
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